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USE OF NUCLEOTIDE ANALOGS IN THE ANALYSIS OF 
OLIGONUCLEOTIDE MIXTURES AND IN HIGHLY MULTIPLEXED NUCLEIC 

ACID SEQUENCING 

Subject matter described herein was developed under NSF Grant 
No. Ger-9452651 . The Government can have certain rights therein. 
RELATED APPLICATIONS 

For U.S. purposes for priority is claimed under 35 U.S.C. §1 19(e) 
to U.S. provisional application Serial No. 60/21 1,356, filed June 13, 
2000, to Charles R. Cantor and Fouad A. Siddiqi, entitled USE OF 
NUCLEOTIDE ANALOGS IN THE ANALYSIS OF OLIGONUCLEOTIDE 
MIXTURES AND IN HIGHLY MULTIPLEXED NUCLEIC ACID 
SEQUENCING." For international purposes benefit of priority is claimed 
thereto. The subject matter of U.S. provisional application Serial No. 
60/21 1,356 is incorporated by reference in it entirety. 

FIELD OF THE INVENTION 

This invention relates to methods, particularly mass spectrometric 
methods, for the analysis and sequencing of nucleic acid molecules. 
DESCRIPTION OF THE BACKGROUND 

Since the recognition of nucleic acid as the carrier of the genetic 
code, a great deal of interest has centered around determining the 
sequence of that code in the many forms in which it occurs. Two studies 
made the process of nucleic acid sequencing, at least with DNA, a 
common and relatively rapid procedure practiced in most laboratories. 
The first describes a process whereby terminally labeled DNA molecules 
are chemically cleaved at single base repetitions (A.M. Maxam and W. 
Gilbert, Proc. Natl. Acad. Sci. USA 74:560-64, 1977). Each base 
position in the nucleic acid sequence is then determined from the 
molecular weights of fragments produced by partial cleavage. Individual 
reactions were devised to cleave preferentially at guanine, at adenine, at 
cytosine and thymine, and at cytosine alone. When the products of these 
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four reactions are resolved by molecular weight, using, for example, 
polyacrylamide gel electrophoresis, DNA sequences can be read from the 
pattern of fragments on the resolved gel. 

In another method DNA is sequenced using a variation of the plus- 
5 minus method (Sanger et al. (1977) Proc. Natl. Acad. Sci. USA 74:5463- 
67, 1977). This procedure takes advantage of the chain terminating 
ability of dideoxynucleoside triphosphates (ddNTPs) and the ability of 
DNA polymerase to incorporate ddNTPs with nearly equal fidelity as the 
natural substrate of DNA polymerase, deoxynucleoside triphosphates 

10 (dNTPs). Briefly, a primer, usually an oligonucleotide, and a template 
DNA are incubated in the presence of a useful concentration of all four 
dNTPs plus a limited amount of a single ddNTP. The DNA polymerase 
occasionally incorporates a dideoxynucleotide that terminates chain 
extension. Because the dideoxynucleotide has no 3'-hydroxyl, the 

15 initiation point for the polymerase enzyme is lost. Polymerization 

produces a mixture of fragments of varied sizes, all having identical 3' 
termini. Fractionation of the mixture by, for example, polyacrylamide gel 
electrophoresis, produces a pattern that indicates the presence and 
position of each base in the nucleic acid. Reactions with each of the four 

20 ddNTPs permits the nucleic acid sequence to be read from a resolved gel. 

These procedures are cumbersome and are limited to sequencing 
DNA. In addition, with conventional procedures, individual sequences are 
separated by, for example, electrophoresis using capillary or slab gels, 
which slow. Mass spectrometry has been adapted and used for 

25 sequencing and detection of nucleic acid molecules (see, e.g., U.S. Patent 
Nos. (6,194,144; 6,225,450; 5,691,141; 5,547,835; 6,238,871; 
5,605,798; 6,043,031; 6,197,498; 6,235,478; 6,221,601; 6,221,605). 
In particular, Matrix-Assisted Laser Desorption/lonization (MALDI) and 
ElectroSpray Ionization (ESI), which allow intact ionization, detection and 
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exact mass determination of large molecules, i.e. well exceeding 300 kDa 
in mass have been used for sequencing of nucleic acid molecules. 

A further refinement in mass spectrometric analysis of high 
molecular weight molecules was the development of time of flight mass 
spectrometry (TOF-MS) with matrix-assisted laser desorption ionization 
(MALDI). This process involves placing the sample into a matrix that 
contains molecules that assist in the desorption process by absorbing 
energy at the frequency used to desorb the sample. Time of flight 
analysis uses the travel time or flight time of the various ionic species as 
an accurate indicator of molecular mass. Due to its speed and high 
resolution, time-of-flight mass spectrometry is well-suited to the task of 
short-range, i.e., less than 30 base sequencing of nucleic acids. Since 
each of the four naturally occurring nucleotide bases dC, dT, dA and dG, 
also referred to herein as C, T, A and G, in DNA has a different molecular 
weight, 

Mc = 289.2 
Mt = 304.2 
Ma = 313.2 
Mg = 329.2, 

where Mc, Mt, Ma, Mg are average molecular weights in daltons of the 
nucleotide bases deoxycytidine, thymidine, deoxyadenosine, and 
deoxyguanosine, respectively, it is possible to read an entire sequence in 
a single mass spectrum. If a single spectrum is used to analyze the 
products of a conventional Sanger sequencing reaction, where chain 
termination is achieved at every base position by the incorporation of 
dideoxynucleotides, a base sequence can be determined by calculation of 
the mass differences between adjacent peaks. In addition, the method 
can be used to determine the masses, lengths and base compositions of 
mixtures of oligonucleotides and to detect target oligonucleotides based 
upon molecular weight. 
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MALDI-TOF mass spectrometry for sequencing DNA using mass 
modification (see, e.g., U.S. Patent Nos. 5,547,835, 6,194,144; 
6,225,450; 5,691,141 and 6,238,871) to increase mass resolution is 
available. The methods employ conventional Sanger sequencing reactions 
5 with each of the four dideoxynucleotides. In addition, for example for 
multiplexing, two of the four natural bases are replaced; dG is substituted 
with 7-deaza-dG and dA with 7-deaza-dA. 

A variety of techniques and combinations thereof have been 
directed to improving the level of accuracy in determining the nucleotide 

10 compositions of mixtures of oligonucleotides using mass spectrometry, 
and many of these methods employ nucleotide analogs. For example, 
Muddiman eta/. (Anal. Chem., 69(8): 1543-1549, 1997) discusses an 
algorithm for the unique definition of the base composition of PCR- 
amplified products, especially longer (> 100bp) oligonucleotides. The 

15 algorithm places a constraint on the otherwise large number of possible 
base compositions for long oligonucleotides by taking into account only 
those masses (measured by electrospray ionization mass spectrometry) 
that are consistent with that of their denatured complementary strands, 
assuming Watson-Crick base-pairing. In addition, the algorithm imposes 

20 the constraint of known primer compositions, since the primer sequences 
are known, and this constraint becomes especially significant with shorter 
PCR products whose mass of "unknown" sequence relative to that of the 
primer mass is small. Muddiman et aL also discusses invoking additional 
measurements for defining the base composition with even greater 

25 accuracy. These include the possibility of post-modifying the PCR 

product using e.g., dimethyl sulfate to selectively methylate every "G" in 
the PCR product, or using a modified base during PCR amplification, 
conducting mass measurements on the modified oligonucleotides, and 
comparing the mass measurements with those of the unmodified 

30 complementary strands. 
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Chen eta/. (Anal. Chem., 71 (1 5): 31 1 8-31 25, 1999) reports a 
method that combines stable isotope 13 C/ 15 N labelling of PCR products 
with analysis of the mass shifts by MALDI-TOF mass spectrometry. The 
mass shift due to labelling of a single type of nucleotide (i.e, A 7 T, G or C) 
5 reveals the number of that type of nucleotide in a given fragment. While 
the method is useful in the measurement and comparison of nucleotide 
compositions of homologous sequences for sequence validation and in 
scoring polymorphisms, tedious repetitive sequencing reactions (using the 
four different labelled nucleotides) and mass spectrometric measurements 
10 are required. 

Hence there is a need in the art for methods that (i) unambiguously 
assign nucleotides in a sequence, and, (ii) resolve large numbers of 
oligonucleotides that have the same length, different base compositions, 
and nearly equal (i.e., less than or equal to about 1 dalton difference) 
15 molecular weights. Therefore it is an object herein to provide methods 
that solve such problems 
SUMMARY OF THE INVENTION 

Provided herein are methods for sequencing and detecting nucleic 
acids using techniques, such as mass spectrometry and gel 
20 electrophoresis, that are based upon molecular mass. The methods use 
deoxynucleotide analogs, modified nucleotide terminators and/or mass- 
labeled primers in one or more reactions for sequencing or detection 
protocols that involve primer extension, and analyze these products from 
more than one oligonucleotide on, for example, a single mass spectrum. 
25 This provides a means for accurate detection and/or sequencing of a an 
oligonucleotide and is particularly advantageous for detecting or 
sequencing a plurality target nucleic acid molecules in a single reaction 
using any technique that distinguishes products based upon molecular 
weight. The methods herein are particularly adapted for mass 
30 spectrometric analyses. 
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For example, a sequencing method provided herein uses 
deoxynucleotide analogs, modified nucleotide terminators and/or mass- 
labeled primers in one or more Sanger sequencing reactions, and analyzes 
these products from more than one oligonucleotide on a single mass 
5 spectrum. In particular, a plurality of primers can be used to 

simultaneously sequence a plurality of nucleic acid molecules or portions 
of the same molecule. By extending the primers with mass-matched 
nucleotides, the resulting products mass shifts that are periodically related 
to the size of the original primer. 

As a result, the sequence of any given oligonucleotide can be 
determined with a high level of accuracy, and also mixtures of a number 
of sequences can be multiplexed in a single mass spectrum. The limit on 
the number of oligonucleotides that can be sequenced simultaneously is 
governed by the base periodicity, the maximum mass shift, and the 
resolving power of analytical tool, such as the mass spectrometer. The 
base periodicity and maximum mass shift can be carefully engineered for 
optimal resolution and accuracy, depending on the number of sequences 
to be simultaneously analyzed, and the information desired; as many 
sequences as desired can be sequenced simultaneously especially in the 
detection and scoring of single nucleotide polymorphisms, insertions, 
deletions and other mutations. 

In another embodiment, a target nucleic acid molecule is 
sequenced using mass-matched nucleotides and chain terminating 
nucleotides. For example, a primer is annealed to a target nucleic acid, 
the primer is extended in the presence of chain-terminating nucleotides 
and mass-matched nucleotides to produce extension products, the 
masses of the extension products follow a periodic distribution that is 
determined by the mass of the mass-matched nucleotides, and the 
sequence of the target nucleic acid is determined from the mass shift of 
each extension product from its corresponding periodic reference mass by 
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virtue of incorporation of the chain terminator. The mass-matched 
nucleotides all have identical masses, and each chain terminating 
nucleotide has a distinct mass that differs from that of the other chain 
terminating nucleotides. This results in unique predetermined values of 
mass shift corresponding to each chain terminating nucleotide and based 
upon the original primer. 

This method is adaptable for any sequencing method or detection 
method that relies upon or includes chain extension. These methods 
include, but are not limited to, sequencing methods based upon Sanger 
sequencing, and detection methods, such as primer oligo base extension 
(PROBE) (see, e.g., U.S. application Serial No. 6,043,031; allowed U.S. 
application Serial No. 09/287,679; and 6,235,478), that rely include a 
step of chain extension. 

Also, contemplated are methods, such as haplotyping methods, in which 
two mutations in the same gene are detection are provided. A detector 
(primer) oligonucleotide is to the hybridized to the first mutation and the 
primer is extended with mass-matched nucleotides and appropriately 
selected chain terminator(s) to detect the second mutation. 

In other embodiments, a plurality of target nucleic acids can be 
multiplexed in a single reaction measurement by annealing each target 
nucl eic acid to a primer of distinct molecular weight each primer is then 
extended with mass-matched nucleotides and chain terminators in 
formats that depend upon whether detection or sequencing is desired. 
These methods are particularly useful for methods of detection in which a 
primer is hybridized to a plurality of target nucleic acid molecules, such as 
immobilized nucleic acid molecules, hybrids separated from unhybridized 
nucleic acids and the detectors detected. Such methods include PROBE, 
in which case the extension reaction is performed in the presence chain 
terminators and mass matched deoxynucleotides. 
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The primers of distinct molecular weight can be selected to differ in 
molecular weight by a value that is greater than the maximum mass shift, 
i.e., the difference in molecular weight between the heaviest and the 
lightest nucleotide terminators in chain extension reactions. The 
5 difference in molecular weight between the primers for a plurality of 
target nucleic acids can be selected to be least 20 daltons greater than 
the maximum mass shift to account for the finite band width of the 
peaks. 

The number of molecules that can be multiplexed is governed by 

10 the periodicity, the maximum mass shift, and the resolving power of the 
sequence detection instrument. In some embodiments, about 7 to about 
25 or more molecules can be multiplexed. For scoring single nucleotide 
polymorphisms, only a single nucleotide terminator is required (depending 
on the base identity of the single nucleotide polymorphism). In this case, 

15 the maximum mass shift required is identically zero, so that larger 
numbers of molecules, greater than 25, 35, 50 and more, can be 
multiplexed, depending on the resolving power of the sequencing format, 
and for mass spectrometry the instrument. Depending on the amount of 
sequence information desired, one, two or three rather than four types of 

20 nucleotide terminators (corresponding to each of the four nucleic acid 
bases) can be used. 

In other embodiments, the mass shift is obtained using pair- 
matched nucleotides, i.e., the mass of each nucleotide base-pair is 
selected so that the masses of all pairs are identical. In one embodiment 

25 thereof, the following steps are performed: (i) the target nucleic acid is 
copied or amplified by a method such as PCR in the presence of the pair- 
matched nucleotide set prior to the sequencing or detection reaction; 
(ii) the target nucleic acid is denatured, and a partially duplex hairpin 
primer is annealed and ligated to the single-stranded template; (iii) the 

30 primer is extended in the presence of chain terminating nucleotides and 
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pair-matched nucleotides to produce extension products, where the 
masses of the extension products follow a periodic distribution that is 
determined by the mass of the pair-matched nucleotide set, and, (iv) the 
target nucleic acid is detected by virtue of its molecular weight or its 
5 sequence is determined from the mass shift of each extension product 
from its corresponding periodic reference mass. 

In another embodiment, the mass of each terminating base pair is 
unique and resolvable, so that the mass shifts corresponding to each 
terminating base pair are unique. The nucleotide terminators are 

10 optionally mass-matched or can be of distinct masses as long as distinct 
values of mass shift are obtained for each terminating base pair. 

In another embodiment, the extension products are treated to 
produce blunt-ended double-stranded extension products by methods 
known to those of skill in the art, such as the use of single-strand specific 

15 nucleases. In an aspect of this embodiment, a plurality of target nucleic 
acids can be multiplexed in a single reaction by annealing each target 
nucleic acid to a primer of distinct molecular weight. The primers can be 
selected to differ in molecular weight by a value that is greater than the 
maximum mass shift, i.e., the difference in molecular weight between the 

20 heaviest and the lightest nucleotide terminating base pairs. Since double 
stranded nucleic acid can be analyzed, the effective sequence read is 
halved relative to the embodiment employing mass-matched nucleotides, 
but the number of molecules that can be multiplexed is doubled, due to 
the increase in period (the value of the mass of a base pair, rather than a 

25 single mass-matched nucleotide). In exemplary embodiments, about 14 
to about 50 sequences are multiplexed. In detection embodiments, about 
50 or more molecules can be simultaneously multiplexed since only a 
single terminating base pair is added in the extension reaction. 

In another embodiment, the chain termination reactions are carried 

30 out separately using a standard nucleotide terminator, pair-matched 
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nucleotides, and mass-labeled primers, if modified nucleotide terminators 
which are either mass-matched or provide distinct values of mass shift for 
each terminating base pair are not available. The reactions are pooled 
prior to detection or sequence analysis. In one embodiment, the mass- 
5 labeled primers can have distinct values of molecular weight that give rise 
to unique values of mass shift or positional mass difference for each 
terminating base. 

In andother method provided herein, a population of nucleic acids 
having the same length but different base compositions can be resolved 

10 by synthesizing the nucleic acids in the presence of a nucleotide analog to 
produce synthesized nucleic acids having incorporated the nucleotide 
analog, where the nucleotide analog is selected to optimally separate the 
masses of the population of nucleic acids according to their individual 
base compositions. For example, the nucleotide analog or analogs are 

15 selected to separate the population of nucleic acids according to base 
composition by greater than 1 dalton. In another embodiment, the 
nucleotide analog or analogs are selected to separate the population of 
nucleic acids according to base composition by mass values of about 3 
daltons to about 8 daltons, depending on the choice of analog and on the 

20 resolving power of the detection instrument. In other embodiments, the 
nucleotide analog or analogs can be selected to restrict oligonucleotides 
having the same length to have the same mass, i.e., a peak separation of 
zero, regardless of differences in base composition, such as in detection 
methods, where it is desirable to separate populations of oligonucleotides 

25 according to their length. 

Nucleic acid molecules that contain mass-matched nucleotides 
and/or pair-matched nucleotides are provided. 

Also provided are combinations for practicing the methods provided 
herein. For instance, in one embodiment, the combinations include a set 

30 of mass-matched deoxynucleotides. In another embodiment, the 
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combinations a set of pair-matched nucleotides and a set of mass- 
matched chain terminating nucleotides. In another embodiment, the 
combination includes a set of pair-matched nucleotides and chain 
terminating nucleotides which form terminating base pairs of distinctly 
5 different molecular weight. In yet another embodiment, the combination 
includes a set of pair-matched nucleotides and mass-labeled primers. In 
other embodiments, mass-staggered primers can be added to as optional 
components. 

Kits containing the combinations with optional instructions and/or 
10 additional reagents are also provided. The kits contain the reagents as 

described herein and optionally any other reagents required to perform the 
reactions. Such reagents and compositions are packaged in standard 
packaging known to those of skill in the art. Additional vials, containers, 
pipets, syringes and other products for sequencing can also be included. 
15 Instructions for performing the reactions can be included. 

Also provided herein are methods for optimization of the analysis of 
base compositions of mixtures of oligonucleotides by mass spectrometry. 
A single spectrum can be used to resolve a very large number of 
oligonucleotides having the same length but different molecular weights 
20 by incorporating a nucleotide analog into the oligonucleotides in the 
mixture such that the peaks are no closer than a minimum value called 
peak separation. The peak separation can be tailored by careful selection 
of the nucleotide analog and of a mass spectrometer with the desired 
resolving power. 

25 The methods herein permit unambiguous and accurate analysis of 

the sequences or molecular weights of large numbers of oligonucleotides 
in a single mass spectrum by combining the rapidity of mass spectrometry 
with the resolving power of nucleotide analogs which are carefully 
selected and incorporated into the oligonucleotide mixture according to 

30 the desired application. 
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Other features and advantages will be apparent from the following 
detailed description and claims. 
BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 shows that when a single spectrum is used to analyze the 
5 products of a conventional Sanger sequencing reaction, where chain 
termination is achieved at every base position by the incorporation of 
dideoxynucleotides, the base sequence can be determined by calculation 
of the mass differences between adjacent peaks (Figures 1a and 1b). 
Figure 2 shows implementation of forced mass modulation using 

10 mass-matched deoxynucleotides. Figure 2a is a simulated mass spectrum 
showing the products and molecular masses of a reaction carried out with 
a suitable polymerase in the presence of a mass-matched nucleotide set 
{"dN") and the four standard dideoxynucleotide terminators. The base 
periodicity is the mass of dN, or 310 daltons. Figure 2b shows a target 

15 second sequence resolved on the same mass spectrum shown in Figure 
2a, using a primer heavier by 77 daltons. The peaks corresponding to the 
reaction products from the first target sequence can fall within the the 
spectrum in Figure 2b, which can never intersect peaks from the second 
target sequence. This permits unambiguous resolution of both sequences 

20 each peak can be uniquely assigned to a nucleotide, a base position, and 
a target sequence. 

Figure 3 shows four different sequences resolved in a single 
spectrum using a set of mass-staggered primers that are separated in 
mass by integer multiples of 77 daltons (77, 154, and 231 daltons). 

25 Figure 4 shows the general implementation of a forced mass 

modulation method using pair-matched nucleotides, for the analysis of 
sequencing reaction products as double-stranded structures. The steps in 
the reaction are as follows: a) a partially duplex hairpin primer with a 3' 
overhang and a 5' phosphate group is annealed and ligated to the single 

30 stranded target sequence; b) the resulting partially duplex structure is 
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subjected to a sequencing reaction using the pair-matched nucleotide set 
described above along with the set of mass-matched terminators (ddM); 
c) products resulting from sequencing reaction b); and, d) the products c) 
from the sequencing reaction are exposed to a strict single strand-specific 
5 nuclease that results in the production of blunt-ended hairpin structures 
ready for analysis by mass spectrometry. 

Figure 5 shows the products and molecular masses of the nuclease 
digestion elucidated in Figure 4d, along with a simulated mass spectrum. 
Figure 6 shows three sequence variants (Figure 6a) that differ from 

10 each other only at a single base position sequenced by a conventional 
Sanger reaction. Figure 6b is a simulated mass spectrum of all reaction 
products shown in Figure 6a. Figure 6c is a graph representing the valid 
sequence permutations that can be elucidated from the mass spectrum 
shown in Figure 6b. Boxed values are fragment masses, solid arrows show 

15 valid sequence branches, dashed arrows represent spurious branches. In 
practice valid branches are indistinguishable from spurious ones. Figure 6d 
is a set of sequences consistent with the graph shown in Figure 6c. 
Spurious sequence reconstructions are shown in lowercase letters, valid 
ones in uppercase letters. 

20 Figure 7a shows the three related sequence variants from Figure 6 

sequenced by Forced Mass Modulation using a single primer and the 
mass-matched nucleotide set from Figure 2 with the standard dideoxy 
terminators. The positions of the differing bases are shown by solid arrows. 
Reaction products are shown along with their respective molecular masses. 

25 Reaction products of variant #2 whose masses differ from those of variant 
#1 are marked with by (*). Reaction products of variant #3 whose masses 
differ from those of variant #1 are marked by ( * *) . Figure 7b is a simulated 
mass spectrum of all reaction products shown in Figure 7a along with a 
sequence graph. The shaded regions represent the only valid mass ranges 

30 that can assumed by the reaction products from Figure 7a. The base 
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periodicity is 310 daltons. Figure 7c is a consensus sequence derived from 
the data shown in Figure 7b, Figure 7d is an expansion of the consensus 
sequence shown in Figure 7c. Spurious reconstructions are shown in 
lowercase letters, valid ones in uppercase letters. Note that there is only a 
5 single spurious reconstruction, as opposed to the eleven errant sequences 
reconstructed from the Sanger reaction described in Figure 6. 

Figure 8 shows the base composition density distributions for the total 
set of possible 7-base oligonucleotides using three different nucleotide sets. 
Note that for the set of naturally occurring bases, nearly every base 

10 composition has its own distinct mass value, but most of these mass values 
are spaced only one dalton from each other. Increasing the peak separation 
also markedly increases the average number of base compositions per 
observed mass, particularly for those masses in the center of the range. 
DETAILED DESCRIPTION OF THE INVENTION 

1 5 Definitions 

Unless defined otherwise, all technical and scientific terms used herein 
have the same meaning as is commonly understood by one of skill in the art 
to which this invention belongs. All patents, patent applications, Genbank 
and other sequence repository sequences, and publications referred to herein 

20 are incorporated by reference. 

As used herein, a biopolymer includes, but is not limited to, nucleic 
acid, proteins, polysaccharides, lipids and other macromolecules. Nucleic 
acids include DNA, RNA, and fragments thereof. Nucleic acids may be 
derived from genomic DNA, RNA, mitochondrial nucleic acid, chloroplast 

25 nucleic acid and other organelles with separate genetic material. 

As used herein "nucleic acid" refers to polynucleotides such as 
deoxyribonucleic acid (DNA) and ribonucleic acid (RNA). The term should 
also be understood to include, as equivalents, derivatives, variants and 
analogs of either RNA or DNA made from nucleotide analogs, single (sense 

30 or antisense) and double-stranded polynucleotides. Deoxyribonucleotides 
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include deoxyadenosine, deoxycytidine, deoxyguanosine and 
deoxythymidine. For RNA, the uracil base is uridine. 

As used herein, "forced mass modulation" refers to methods provided 
herein that use deoxynucleotide analogs, modified nucleotide terminators, 
5 mass-labeled primers, mass-staggered primers and other such nucleotides, 
nucleic acids and analogs thereof, to unambiguously assign peak positions 
of mass fragments of oligonucleotides according to their base position, base 
identity, and target sequence from which the fragments arose. The method 
is used to sequence, detect or identify single oligonucleotide or plurality 

10 thereof. Hence the method is used, for example for muliplex sequencing and 
detection of nucleic acid molecules among mixtures thereof. 

As used herein, "nucleotides" include, but are not limited to, the 
naturally occurring nucleoside mono-, di-, and triphosphates: deoxyadenosine 
mono-, di- and triphosphate; deoxyguanosine mono-, di- and triphosphate; 

15 deoxythymidine mono-, di- and triphosphate; and deoxycytidine mono-, di- 
and triphosphate (referred to herein as dA, dG, dT and dC or A, G, T and C, 
respectively). Nucleotides also include, but are not limited to, modified 
nucleotides and nucleotide analogs such as deazapurine nucleotides, e.g., 
7-deaza-deoxyguanosine (7-deaza-dG) and 7-deaza-deoxyadenosine (7- 

20 deaza-dA) mono-, di- and triphosphates, deutero-deoxythymidine (deutero- 
dT) mon-, di- and triphosphates, methylated nucleotides e.g., 5- 
methyldeoxycytidine triphosphate, 13 C/ 15 N labelled nucleotides and 
deoxyinosine mono-, di- and triphosphate. For those skilled in the art, it will 
be clear that modified nucleotides and nucleotide analogs can be obtained 

25 using a variety of combinations of functionality and attachment positions. 

As used herein, a complete set of chain-elongating nucleotides refers 
to the four different nucleotides or analogs thereof that hybridize to each of 
the four different bases comprising the nucleic acid template. 

As used herein, the term "mass-matched nucleotides" refers to a set 

30 of nucleotide analogs wherein each analog is of identical mass to each of the 




other analogs. For example, analogs of dA, dG, dC and dT can form a mass- 
matched nucleotide set, when each analog is selected to have the same 
molecular weight as the others in the set. Mass-matched nucleotide sets 
can be identified by selecting chemically modified derivatives of natural 
5 bases or by the use of a universal base analog such as deoxyinosine or 5- 
nitroindole and 3-nitropyrrole (5-nitroindole and 3-intropyrrole can be in the 
dideoxy form) which can form base pairs with more than one of the natural 
bases. Others include, 3-methyl 7-propynyl isocarbostyril, 5-methyl 
iscarbostyril, and 3-methyl iscarbostyril. As a result, oligonucleotides that 

10 contain such bases differ in molecular weight only as a function of length 
thereof. Furthermore, incorporation of a single nucleotide(s) that is (are) not 
in the set renders such the oligonucleotide(s) readily identifiable by mass, 
particularly by spectrometric analysis. 

As used herein, the term "pair-matched nucleotides" refers to a 

15 nucleotide set in which the nucleotide analogs are selected such that the 
total mass each base pair is identical. For example, replacing dG with the 
nucleotide analog 7-deaza-dG forces the mass of each base pair, i.e., (dA + 
dT) and (dC + 7-deaza-dG) to be identical. Exemplary pair-matched 
nucleotides, include, but are not limited to, 7-deaza-dA + 

20 phosphorothioate-dT ((31 2.2 + 320.2) = 632.4 Da) and 5-methyl-dC + dG 
((303.2 + 329.2) = 632.4 Da); phosphorothioate-7-deaza-dA + dU 
((328.2 + 290.2) - 61 8.4 Da) and dC + dG = ((289.2 + 329.2) - 618.4 
Da), and other such pairs that may be readily selected. Another exemplary 
set of mass-matched nucleotides with a molecular mass of 328.2: 

25 7-deaza-dG, phosphorothioate-7-deaza-dA, 5-propynyl-dU and 5-cyano- 
methyl-2'-deoxycytidine. 

As used herein, the term "nucleotide terminator" or "chain terminating 
nucleotide" refers to a nucleotide analog that terminates nucleic acid 
polymer (chain) extension during procedures wherein a DNA template is 

30 being sequenced or replicated. The standard chain terminating nucleotides, 
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i.e., nucleotide terminators include 2',3 / -dideoxynucleotides(dd ATP, ddGTP, 
ddCTP and ddTTP, also referred to herein as dideoxynucleotide terminators) . 
As used herein, dideoxynucleotide terminators also include analogs of the 
standard dideoxynucleotide terminators, e.g. , 5-bromo-dideoxyuridine, 5- 
5 methyl-dideoxycytidine and dideoxyinosine are analogs of ddTTP, ddCTP and 
ddGTP, respectively. 

As used herein, "mass-matched terminators" refers to a set of 
nucleotide terminators that are selected such that each analog of ddA, ddG, 
ddC and ddT making up the mass-matched set has exactly the same 

10 molecular weight. Mass-matched terminator sets can be constructed by 
selecting chemically modified derivatives of standard dideoxynucleotides or 
by the use of a universal dideoxynucleotide analog that form base pairs with 
more than one of the natural bases. Exemplary mass-matched nucleotides 
include, but are not limited to, 3-methyl 7-propynyl isocarbostyril, 5-methyl 

15 iscarbostyril and 3-methyl iscarbostyril. As used herein, the terms 
"oligonucleotide" or "nucleic acid" refer to single-stranded and/or double- 
stranded polynucleotides such as deoxyribonucleic acid (DNA), and 
ribonucleic acid (RNA) as well as derivatives of either RNA or DNA into 
which nucleotide or dideoxynucleotide analogs have been incorporated. Also 

20 included in the term "nucleic acid" are analogs of nucleic acids such as 
peptide nucleic acid (PNA), phosphorothioate DNA, and other such analogs 
and derivatives. 

As used herein, "nucleotide composition" or "base composition" 
refers to the numerical ratio of the four nucleotide bases relative to each 
25 other in an oligonucleotide. 

As used herein, a target nucleic acid refers to any nucleic acid of 
interest in a sample. It can contain one or more nucleotides. A target 
nucleotide sequence refers to a particular sequence of nucleotides in a target 
nucleic acid molecule. Detection or identification of such sequence results 
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in detection of the target and can indicate the presence or absence of a 
particular mutation or polymorphism. 

As used herein, "partially duplex hairpin" refers to a partially self- 
complementary oligonucleotide, which forms intramolecular base-pairs within 
5 its self-complementary region, leaving a "loop" of bases at one end of the 
molecule and a single-stranded "overhang" region at the other end. Thus, 
the oligonucleotide assumes a hairpin-like motif. "Blunt-ended hairpin 
structures", as referred to herein, are similar to the partially duplex hairpin 
structures with the exception that they do not have a single-stranded 

10 "overhang" region. 

As used herein, "base periodicity" or "period" (Pbase) refers to the 
quasi-periodic distribution of the molecular weights of products obtained 
using Forced Mass Modulation. The base periodicity results from either the 
mass of the mass-matched deoxynucleotide set, or the mass of the pair- 

15 matched deoxynucleotide set, or from the modified chain terminators 
depending on the embodiment implemented. The base sequence or nucleic 
acid molecule identity is encoded in the pattern (or detectable therein) in 
which the observed mass distribution deviates from absolute regular 
periodicity. 

20 As used herein, the "periodic reference mass" at base position "n" 

in any given oligonucleotide molecule, MpR[n], is defined as the sum of: (i) 
the mass of the primer (Mprimer) used to sequence the DNA template using 
Forced Mass Modulation, (ii) the mass of the lightest nucleotide terminator 
(Mught), and, (iii) (n-1) multiple of the base periodicity Pbase. 

25 As used herein, the "positional mass difference" or "mass shift" at 

base position "n" in any given oligonucleotide molecule, Mdiff[n], is defined 
as the distance in daltons between the observed peak, Mobsfn], and the nth 
periodic reference mass. 

As used herein, the "maximum mass shift" Smax is the maximum 

30 possible value of the positional mass difference, depending on the choice of 
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mass-matched nucleotides and nucleotide terminators used in the 
implementation of Forced Mass Modulation. Accordingly, the maximum 
mass shift can be modulated by the choice of mass-matched nucleotides and 
nucleotide terminators. 
5 As used herein, a "primer" refers to an oligonucleotide that is 

suitable for hybridizing, chain extension, amplification and sequencing. 
Similarly, a probe is a primer used for hybridization. The primer refers to a 
nucleic acid that is of low enough mass, typically about between about 5 
and 200 nucleotides, generally about 70 nucleotides or less than 70, and of 

10 sufficient size to be conveniently used in the methods of amplification and 
methods of detection and sequencing provided herein. These primers 
include, but are not limited to, primers for detection and sequencing of 
nucleic acids, which require a sufficient number nucleotides to form a stable 
duplex, typically about 6-30 nucleotides, about 10-25 nucleotides and/or 

15 about 1 2-20 nucleotides. Thus, for purposes herein, a primer is a sequence 
of nucleotides contains of any suitable length, typically containing about 6- 
70 nucleotides, 1 2-70 nucleotides or greater than about 1 4 to an upper limit 
of about 70 nucleotides, depending upon sequence and application of the 
primer. 

20 As used herein, the term "mass-labeled primers" refers to a set of 

primers that differ in mass by values that provide distinct and resolvable 
positional mass differences for each of the four termination reactions in an 
embodiment of Forced Mass Modulation. In this particular embodiment of 
Forced Mass Modulation, each of the termination reactions for a given 

25 oligonucleotide is carried out separately using each of the mass-labeled 
primers, and the reaction products are combined prior to obtaining a mass 
spectrum. 

As used herein, the term "mass-staggered primers" refers to the mass 
difference ("staggering" of the masses) between the primers used in 
30 multiplexed sequencing using Forced Mass Modulation. For resolution of 
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multiple sequences using this method, the differences between the masses 
of the primers should at least be equal to the maximum mass shift, and is 
generally greater than the maximum mass shift by at least 20 daltons to 
account for the finite width of each observed peak. 
5 As used herein, reference to mass spectrometry encompasses any 

suitable mass spectrometric format known to those of skill in the art. Such 
formats include, but are not limited to, Matrix-Assisted Laser 
Desorption/lonization, Time-of-Flight (MALDI-TOF), Electrospray (ES), IR- 
MALDI (see, e.g., published International PCT application No. 99/5731 8 and 

10 U.S. Patent No. 5,118,937), Ion Cyclotron Resonance (ICR), Fourier 
Transform and combinations thereof. MALDI, particular UV and IR, are 
among the preferred formats. 

As used herein, mass spectrum refers to the presentation of data 
obtained from analyzing a biopolymer or fragment thereof by mass 

15 spectrometry either graphically or encoded numerically. 

As used herein, used herein, pattern with reference to a mass 
spectrum or mass spectrometric analyses, refers to a characteristic 
distribution and number of signals (such peaks or digital representations 
thereof). 

20 As used herein, signal in the context of a mass spectrum and analysis 

thereof refers to the output data, which the number or relative number of 
molecules having a particular mass. Signals include "peaks" and digital 
representations thereof. 

As used herein, "mass spectrum division multiplexing" is an 

25 embodiment of Forced Mass Modulation in which unambiguous resolution 
of multiple sequences in a single spectrum is possible by judicious selection 
of mass staggered primers. 

As used herein, "analysis" refers to the determination of certain 
properties of a single oligonucleotide, or of mixtures of oligonucleotides. 

30 These properties include, but are not limited to, the nucleotide composition 
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and complete sequence of an oligonucleotide or of mixtures of 
oligonucleotides, the existence of single nucleotide polymorphisms between 
more than one oligonucleotide, the masses and the lengths of 
oligonucleotides and the presence of a molecule or sequence within 
5 molecule in a sample. 

As used herein, "multiplexing" refers to the simultaneous 
determination of more than one oligonucleotide molecule, or the 
simultaneous analysis of more than one oligonucleotide, in a single mass 
spectrometric or other sequence measurement, i.e., a single mass spectrum 
10 or other method of reading sequence. 

As used herein, "polymorphisms" refer to variants of a gene or an 
oligonucleotide molecule that differ at more than one base position. In 
"single nucleotide polymorphisms", the variants differ at only a single base 
position. 

15 As used herein, amplifying refers to means for increasing the amount 

of a bipolymer, especially nucleic acids. Based on the 5' and 3' primers that 
are chosen, amplification also serves to restrict and define the region of the 
genome which is subject to analysis. Amplification can be by any means 
known to those skilled in the art, including use of the polymerase chain 

20 reaction (PCR) etc. Amplification, e.g., PCR must be done quantitatively 
when the frequency of polymorphism is required to be determined. 

As used herein, "polymorphism" refers to the coexistence of more 
than one form of a gene or portion thereof. A portion of a gene of which 
there are at least two different forms, i.e., two different nucleotide 

25 sequences, is referred to as a "polymorphic region of a gene". A 
polymorphic region can be a single nucleotide, the identity of which differs 
in different alleles. A polymorphic region can also be several nucleotides in 
length. Thus, a polymorphism, e.g. genetic variation, refers to a variation 
in the sequence of a gene in the genome amongst a population, such as 

30 allelic variations and other variations that arise or are observed. Thus, a 
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polymorphism refers to the occurrence of two or more genetically 
determined alternative sequences or alleles in a population. These 
differences can occur in coding and non-coding portions of the genome, and 
can be manifested or detected as differences in nucleic acid sequences, gene 
5 expression, including, for example transcription, processing, translation, 
transport, protein processing, trafficking, DNA synthesis, expressed proteins, 
other gene products or products of biochemical pathways or in post- 
translational modifications and any other differences manifested amongst 
members of a population. A single nucleotide polymorphism (SNP) refers to 

10 a polymorphism that arises as the result of a single base change, such as an 
insertion, deletion or change in a base. 

A polymorphic marker or site is the locus at which divergence occurs. 
Such site may be as small as one base pair (an SNP). Polymorphic markers 
include, but are not limited to, restriction fragment length polymorphisms, 

15 variable number of tandem repeats (VNTR's), hypervariable regions, 
minisatellites, dinucleotide repeats, trinucleotide repeats, tetranucleotide 
repeats and other repeating patterns, simple sequence repeats and 
insertional elements, such as Alu. Polymorphic forms also are manifested 
as different mendelian alleles for a gene. Polymorphisms may be observed 

20 by differences in proteins, protein modifications, RNA expression 
modification, DNA and RNA methylation, regulatory factors that alter gene 
expression and DNA replication, and any other manifestation of alterations 
in genomic nucleic acid or organelle nucleic acids. 

As used herein, "polymorphic gene" refers to a gene having at least 

25 one polymorphic region. 

As used herein, "allele", which is used interchangeably herein with 
"allelic variant" refers to alternative forms of a gene or portions thereof. 
Alleles occupy the same locus or position on homologous chromosomes. 
When a subject has two identical alleles of a gene, the subject is said to be 

30 homozygous for the gene or allele. When a subject has two different alleles 
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of a gene, the subject is said to be heterozygous for the gene. Alleles of a 
specific gene can differ from each other in a single nucleotide, or several 
nucleotides, and can include substitutions, deletions, and insertions of 
nucleotides. An allele of a gene can also be a form of a gene containing a 
5 mutation. 

As used herein, "predominant allele" refers to an allele that is 
represented in the greatest frequency for a given population. The allele or 
alleles that are present in lesser frequency are referred to as allelic variants. 
As used herein, a subject, includes, but is not limited to, animals, 
10 plants, bacteria, viruses, parasites and any other organism or entity that has 
nucleic acid. Among subjects are mammals, preferably, although not 
necessarily, humans. A patient refers to a subject afflicted with a disease 
or disorder. 

As used herein, a phenotype refers to a set of parameters that 
15 includes any distinguishable trait of an organism. A phenotype can be 
physical traits and can be, in instances in which the subject is an animal, a 
mental trait, such as emotional traits. 

As used herein, "resolving power" of a mass spectrometer is the ion 
separation power of the instrument, i.e., it is a measure of the ability of the 
20 mass spectrometer to separate peaks representing different masses. The 
resolving power R is defined as m/Am, where m is the ion mass and Am is 
the difference in mass between two resolvable peaks in a mass spectrum. 

As used herein, "assignment" refers to a determination that the 
position of a nucleic acid fragment indicates a particular molecular weight 
25 and a particular terminal nucleotide. 

As used herein, "a" refers to one or more. 

As used herein, "plurality" refers to two or more, up to an amount 
that is governed by the base periodicity, the maximum mass shift, and the 
resolving power of the mass spectrometer. 
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As used herein, an array refers to a pattern produced by three or more 
items, such as three or more loci on a solid support. 

As used herein, "distinct" refers to a unique value of molecular 
weight, mass shift or period that is different from every other value of 
5 molecular weight, mass shift or period in the measurement. 

As used herein, "unambiguous" refers to the unique assignment of a 
particular oligonucleotide fragment according to the identity of its terminal 
base position and, in the event that a number of molecules are multiplexed, 
that the peak representing an oligonucleotide fragment can also be uniquely 

10 assigned to a particular molecule. 

As used herein, the symbols Mc, Mt, Ma and Mg are average 
molecular weights in daltons of the nucleotides deoxycytidine, thymidine, 
deoxyadenosine and deoxyguanosine, respectively, or of analogs thereof. 
Mavg, the average molecular weight of any given oligonucleotide is a function 

15 of the average molecular weights of each of the nucleotides comprising the 
oligonucleotide, the numbers c, t, a and g of each nucleotide present in the 
oligonucleotide, the length of the oligonucleotide n' that is the sum of c, t, 
a and g, and the constant k that represents the mass of any other chemical 
groups on the molecule, such as terminal phosphates. 

20 As used herein, Ntotal is the total number of possible base 

compositions for an oligonucleotide of length n\ 

As used herein, "peak separation" or "minimum peak separation" S 
refers to the minimum value of the distance between consecutive peaks in 
a mass spectrum that resolves a large number of oligonucleotides having the 

25 same lengths but different molecular weights, i.e., different base 
compositions. The peak separation, which can be tailored by careful 
selection of the nucleotide analogs incorporated into the oligonucleotide and 
by a mass spectrometer of desired resolving power, is usually a positive 
integer greater than one, and typically a positive integer greater than or 

30 equal to 3. For two oligonucleotides having the same length n' but different 
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base compositions, their molecular weights will either correspond to the 
same peak if the molecular weights are identical, or to two peaks separated 
at least by a value equal to the peak separation. 

As used herein, L is the maximum number of allowed oligonucleotide 
5 masses for a given nucleotide set. It is directly proportional to the 
oligonucleotide length n' and the mass difference between the heaviest and 
lightest nucleotides in the set, and is inversely proportional to the peak 
separation. 

As used herein, D refers to the average density of different base 
10 compositions per allowed mass value, given the set of all possible base 
compositions of an oligonucleotide of length n\ 

As used herein, Mheavy refers to the mass of the heaviest nucleotide, 
nucleotide terminator or terminating base pair in daltons, depending on the 
specific embodiment of Forced Mass Modulation being described. 
15 As used herein, Miight refers to the mass of the lightest nucleotide, 

nucleotide terminator or terminating base pair in daltons, depending on the 
specific embodiment of Forced Mass Modulation being described. 
As used herein, Mphmer is the mass of the primer in daltons. 
As used herein, Mobs[n] is the observed mass of the sequencing 
20 reaction at the nth base position. 

As used herein, Mterm[n] refers to the mass in daltons of the nth 
terminating nucleotide. 

As used herein, L' is the theoretical upper limit on the number of 
sequences that be multiplexed in a single mass spectrum. L' is directly 
25 proportional to the base periodicity Pbase, and is inversely proportional to the 
maximum mass shift Smax. 

As used herein, Mdupiex is the mass in daltons of the fully duplex hairpin 
primer in the implementation of Forced Mass Modulation using pair-matched 
nucleotides. 
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As used herein, MddM is the mass in daltons of a dideoxy terminator 
that belongs to a set of mass-matched terminators. 

As used herein, Mtar 9 [n] is the mass of the nth nucleotide past the 
priming site in the 3' to 5' direction in the target sequence, i.e., the 
5 oligonucleotide whose sequence is being determined. 

As used herein, "specifically hybridizes" refers to hybridization of a 
probe or primer only to a target sequence preferentially to a non-target 
sequence. Those of skill in the art are familiar with parameters that affect 
hybridization; such as temperature, probe or primer length and composition, 
10 buffer composition and salt concentration and can readily adjust these 
parameters to achieve specific hybridization of a nucleic acid to a target 
sequence. 

As used herein, a biological sample refers to a sample of material 
obtained from or derived from biological material, such as, but are not limited 
15 to, body fluids, such blood, urine, cerebral spinal fluid and synovial fluid, 
tissues and organs. Derived from means that sample can be processed, 
such as by purification or isolation and/or amplification of nucleic acid 
molecules. 

As used herein, a composition refers to any mixture. It may be a 
20 solution, a suspension, liquid, powder, a paste, aqueous, non-aqueous or 
any combination thereof. 

As used herein, a combination refers to any association between two 

or among more items. 

As used herein, "kit" refers to a package that contains a combination 
25 and optionally instructions and/or reagents and apparatus for use with the 
combination. 
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Forced Mass Modulation for analysis of nucleic acid molecules 
Time of flight analysis and drawbacks thereof 

While time-of-flight mass spectrometry offers a number of advantages 
over conventional techniques such as gel electrophoresis, the peculiar 
relationship between the masses of the bases in DNA complicates the 
analysis of complex mixtures of oligonucleotides by mass spectrometry. For 
a given oligonucleotide, the average molecular weight, Mavg, is given by the 
following equation: 

(i) Mavg = k + cMc + tMT + aMA + gMc 

where Mc, Mt, Ma, Mg are the average molecular weights of each of the four 
nucleotide bases (cytosine, thymine, adenine, guanine) and c, t, a, g 
represent the number of each base present in the oligonucleotide. The term 
k is a constant representing the mass of any other chemical groups on the 
molecule, such as terminal phosphates. Rearranging equation (i) to give the 
average molecular weight as a function of the length of the oligonucleotide 
in bases yields 

(ii) Mavg = k + n'Mc + KMt - Mc) + a(MA - Mc) + g(Ms - Mc) 
where n', the oligonucleotide length, is defined as 

n'=c+t + a + g 
Substituting the masses of the naturally occurring bases in DNA (to one- 
tenth dalton): 

Mc = 289.2 
Mt = 304.2 
Ma = 313.2 
Mg = 329.2 

into equation (ii) yields 

(iii) Mavg = k + 289.2n' + t(304.2 - 289.2) + a(31 3.2 - 289.2) + 

g(329.2 - 289.2), 
which can be simplified to 

(iv) Mavg = k + 289. 2n' + 1 5t + 24a + 40g 
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Close inspection of equation (iv) reveals that it is almost always possible to 
find two oligonucleotides of the same length but of different base 
composition whose average masses differ by only one dalton. For example, 
all 7-mers having a base composition of A2C2G2T have an average molecular 
5 weight of (21 67.4 + k), while all 7-mers with the base composition A3CGT2 
have an average molecular weight of (2166.4 + k). Since the following 
relation 

(Mc + Mg) = (Mt + Ma) + 1 
is always true for the naturally occurring bases in DNA, simply replacing one 

10 C and one G in an oligonucleotide with one A and one T will produce a new 
oligonucleotide exactly one dalton lighter. Many other "single-dalton 
difference" relations, such as 

4Ma = (Mc + Mt + 2Mg) + 1 
can readily be found for the naturally occurring bases. 

1 5 Thus, the possibility always exists that two or more oligonucleotides 

of same length and different molecular weight (and, therefore, different base 
composition) will be too close in mass to be resolved by a time-of-flight 
instrument. Two oligonucleotides of same length but different molecular 
weight differ in base composition unless they are each composed of different 

20 nucleotide analogs, whereas two oligonucleotides of same length and same 

molecular weight can have either the same or different base compositions. 

This problem becomes increasingly severe with increasing oligonucleotide 

size, since the total number of possible base compositions, Ntotal, scales as 

a cubic function of the oligonucleotide length n', in bases: 

25 (v) Ntotal = In' + 1)(n' + 2)(n 7 + 3) 

6 

The use of time-of-flight mass spectrometry in sequencing applications 
also poses several potential problems. The great drawback of sequencing 
by the Sanger method is that the molecular weights of the Sanger reaction 
30 products can appear virtually anywhere on the mass axis depending on the 
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particular sequence being examined. As a result, the absolute mass of any 
single Sanger fragment has to 

be measured with sufficient accuracy to calculate its distance from the 
masses from the fragments above and below it. Thus, determination of the 
5 identity of a single base depends on the accuracy of two separate mass 
measurements. Any error in a determination mass of a single fragment 
affects the accuracy of two bases in the sequence. 

For longer sequences (30-50 bases), it may not be possible to 
determine the mass difference between adjacent peaks with sufficient 

10 accuracy to unambiguously determine base identity. This is particularly a 
problem for the nucleotides A and T 7 which differ in mass only by nine 
daltons. The problem is addressed by resolving each of the four termination 
reactions in a separate mass spectrum. In this case each peak functions 
essentially as binary signal indicating the presence of a base at a particular 

15 position, much as in conventional electrophoretic sequencing. Using 
separate spectra, however, increases read accuracy but at the expense of 
increasing the number of required mass measurements by a factor of four. 

It is possible to resolve two target sequences by the Sanger method 
in a single mass spectrum, provided that all products of the sequencing 

20 reactions have unique and resolvable masses, and multiplex methods using 
mass modified bases have been developed. But, where two or more 
reaction products have the same mass, then unambiguous reconstruction of 
the two target sequences is not possible (see, e.g., Figures 1c-e). In 
addition, there is no way to determine a priori which observed masses 

25 belong to a particular sequence. In practice, this means that multiplexed 
Sanger sequencing by mass spectrometry can be difficult. The methods 
provided herein resolve these problem and provide a way to determine a 
priori which masses are associated with extension of a particular primer. 
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Forced Mass Modulation 

As noted above. Forced Mass Modulation refers to methods provided 
herein that permit unambiguously assign peak positions (or masses) to mass 
fragments of oligonucleotides according to their base position, base identity, 
5 and target sequence from which the fragments arose. The methods use 
deoxynucleotide analogs, modified nucleotide terminators, mass-labeled 
primers, mass-staggered primers and other such nucleotides, nucleic acids 
and analogs thereof to provide a means for deconvoluting complex mass 
spectra or output from other mass determining techniques. These methods 

1 0 permit deconvolution of highly multiplexed nucleic acid reaction mixtures for 
sequencing methods and detection methods that include a step of primer 
extension. In practicing these methods, primers are extended using mass- 
matched nucleotides and chain terminators (or in some embodiments mass 
where it is only necessary to detect incorporation (or the absence of 

1 5 incorporated) mass-matched terminators and optionally mass-matched chain 
extending nucleotides) . Because the sequence and/or molecular mass of a 
primer is known, and the extended nucleotides have the same molecular 
mass, a periodicity in molecular mass that is a function of molecular weight 
of the selected mass matched nucleotide(s) results. 

20 As described in more detail below, for sequencing reactions using 

chain terminators, the deviation from the periodicity results from 
incorporation of a chain terminator. The deviation is a function of the 
particular terminator incorporated. For detection methods, incorporation of 
a terminator will indicate the presence of a mutation (if the terminator is 

25 selected to pair with the first mutated nucleotide. Any shift from periodicity 
will indicate the presence of the mutation. These methods, thus provide a 
simple, reliable way to detect the presence of a mutations or target 
nucleotide(s) in a sequence and to sequence nucleic acids. Forced mass 
modulation can be used with any method, such as mass spectrometry and 
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gel electrophoresis, that relies on molecular weight as an output. Mass 
spectrometry is exemplified herein. 

The methods, designated Forced Mass Modulation methods, provided 
herein, are implemented by suitable selection of nucleotides and/or chain 
5 terminators, such as by the use of deoxynucleotide analogs, modified 
nucleotide terminators and mass-labeled primers in one or more reactions. 
Forced Mass Modulation can be used to simultaneously sequence or detect 
large numbers, such twenty-five or more) oligonucleotides, with a high 
degree of resolution and accuracy. It can also be used to simplify the 

10 analysis of closely related sequence variants, as is required in the detection 
and scoring of nucleotide polymorphisms, including single nucleotide 
polymorphism (SNPs) and for other genotypical analyses. Forced Mass 
Modulation greatly improves the use of mass spectrometry for nucleic acid 
analyses. Nearly every application relies on mass measurements that can 

15 benefit in increased accuracy and in a reduction of the number of required 
spectra. Another advantage of Forced Mass Modulation is the number of 
different ways in which it can be implemented, allowing it to be tailored to 
particular experimental or instrumental limitations. 

For example, compared to the conventional Sanger methods, Forced 

20 Mass Modulation, provides increased accuracy, simplified interpretation of 
mass data, and the ability to use a single mass spectrum for the 
unambiguous resolution of several distinct nucleic acid molecules. For mass 
spectrometry applications, the methods provide unambiguous assignment of 
peak positions of mass fragments of oligonucleotides according to their base 

25 position, base identity, and target sequence from which the fragments arose. 
Thus, the methods herein are advantageously used for multiplexing, in 
which a plurality of reactions are run in a single reaction (single pot). Forced 
Mass Modulation, exemplified with reference to sequencing methods, such 
as PROBE, can also be adapted to detection methods in which a primer is 

30 extended. 
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In Forced Mass Modulation in which a primer is extended with mass- 
matched nucleotides, for examples, the molecular weights of extended 
nucleic acid chains, such as sequencing reaction products, are constrained 
since all extension products from the same primer will have a molecular 
5 weight that differs either by the length of the extension and the chain 
terminator. As a result, the extension products assume a quasi-periodic 
distribution on the mass axis with a predetermined base periodicity. For 
sequencing, the base sequence itself is encoded in the pattern in which the 
observed mass distribution deviates from absolute regular periodicity. Since 

10 the base periodicity will always be known a priori, since the primer is known, 
each peak in the observed mass spectrum can be matched unambiguously 
to a unique nucleotide position in the target sequence. The initiating primers 
fix each set of nested fragments or extended products, and the use of mass- 
matched nucleotides creates the periodicity, 

15 As demonstrated by the Examples below, the method is advantageous 

for numerous applications including sequencing and a variety of detection 
methods, including primer oligo base extension (PROBE) (see, e.g., U.S. 
application Serial No. 6,043,031; allowed U.S. application Serial No. 
09/287,679; and 6,235,478) that use mass spectrometry to distinguish 

20 between extended primers. If the base compositions of the target 
oligonucleotides are known a priori then it is possible to select a nucleotide 
set that produces oligonucleotide masses that are distinct and resolvable for 
any particular instrument or application. 

Conversely, it is also possible to select a nucleotide set that restricts 

25 specific oligonucleotides to have the same mass, regardless of a change in 
base composition. The strategy of restricting specific oligonucleotides to 
have the same mass can be used to separate more than one oligonucleotide 
population of different lengths by restricting all oligonucleotides of a 
particular length to the same molecular weight, irrespective of differences in 

30 base composition. 
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The oligonucleotide analysis or sequencing in methods provided herein 
can be accomplished by one of several methods employed in the art for the 
synthesis, resolution and/or detection of nucleic acids. Depending on the 
embodiment implemented, modified nucleotides can be incorporated into the 
5 oligonucleotides by chemical (Oligonucleotides and Analogues: A Practical 
Approach, F. Eckstein, ed., IRL Press Oxford, 1 991 ) or enzymatic (F. Sanger 
et al., Proc. Natl. Acad. Sci. USA 74:5463-67, 1977) synthesis. Extension 
products or truncated products of the oligonucleotides to be sequenced can 
be obtained using chemical (A.M. Maxam and W. Gilbert, Proc. Natl. Acad. 

10 Sci. USA 74:560-64, 1 977) or enzymatic (F. Sanger et al., Proc. Natl. Acad. 
Sci. USA 74:5463-67, 1977) methods. 

For the resolution and detection of target nucleic acids any mass 
determination method, such as, but are not limited to, chromatography, gel 
electrophoresis, capillary zone electrophoresis and mass spectrometry, is 

1 5 used. Mass spectrometric formats, include, but are not limited to, are matrix 
assisted laser desorption ionization (MALDI), electrospray (ES), ion cyclotron 
resonance (ICR) and Fourier Transform. For ES, the samples, dissolved in 
water or in a volatile buffer, are injected either continuously or 
discontinuously into an atmospheric pressure ionization interface (API) and 

20 then mass analyzed by a quadrupole. The generation of multiple ion peaks 
which can be obtained using ES mass spectrometry can increase the 
accuracy of the mass determination. Even more detailed information on the 
specific structure can be obtained using an MS/MS quadrupole configuration 
in MALDI mass spectrometry, various mass analyzers can be used, 

25 e.g., magnetic sector/magnetic deflection instruments in single or triple 
quadrupole mode (MS/MS), Fourier transform and time-of-flight (TOF) 
configurations as is known in the art of mass spectrometry. For the 
desorption/ionization process, numerous matrix/laser combinations can be 
used. Ion-trap and reflectron configurations can also be employed. 
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Pair-matched nucleotide-based methods 

Forced Mass Modulation can be implemented using a deoxynucleotide 
set in which the mass of each base pair is identical, termed a pair-matched 
nucleotide set. A pair-matched nucleotide set can easily be formed, for 
5 example, by replacing dG (329.2 Da) in the set of naturally occurring 
nucleotides with 7-deaza-dG (328.2 Da). This forces the mass of each base 
pair to be 617.4 daltons: 

(dA + dT) = (313.2 + 304.2) = 617.4 Da 

(dC + 7-deaza-dG) = (289.2 + 328.2) = 617.4 Da 

10 Many other pair-matched sets are possible using available nucleotide 
analogs. For this embodiment, the target DNA sequence can be 
composed entirely of the pair-matched nucleotide set. This can be 
accomplished by amplifying the target DNA sequence by PCR using the 
pair-matched nucleotide set prior to the sequencing reaction. 

15 a further requirement for this embodiment of Forced Mass 

Modulation is that the mass of each terminating base pair is unique and 
resolvable. The standard dideoxy terminators therefore cannot be used 
with the pair-matched nucleotide set described above, because the 
masses of all terminating base pairs are identical at 601 .4 daltons, except 

20 ddG:dC, which is 602.4 daltons. For the sake of clarity in this example, 
it is assumed that a set of mass-matched terminators is available ("ddM," 
defined as set of chain-terminating nucleotides that have exactly the same 
molecular weight ddA = ddC = ddG = ddT). If the mass of ddM is 
arbitrarily chosen to be 500 daltons, then the masses of the terminating 

25 base pairs are as follows: 
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Terminating Base Pair 



Mass (Da) 



ddM: dC 



789.2 



ddM: dT 



804.2 



ddM: dA 



813.2 



5 



ddM: 7-deaza-dG 



828.2 



In practice it is also possible to implement Forced Mass Modulation using 
a set of terminators that have different masses, this is discussed in detail 
below. 

Exemplary embodiments in which the mass shift is obtained using 
10 pair-matched nucleotides, where the mass of each nucleotide base-pair is 
selected so that the masses of all pairs are identical, are described in 
Example 4. In one embodiment thereof, the following steps are 
performed: (i) the target nucleic acid is copied or amplified by a method 
such as PCR in the presence of the pair-matched nucleotide set prior to 
15 the sequencing or detection reaction; (ii) the target nucleic acid is 

denatured, and a partially duplex hairpin primer is annealed and ligated to 
the single-stranded template; (iii) the primer is extended in the presence of 
chain terminating nucleotides and pair-matched nucleotides to produce 
extension products; (iv) the masses of the extension products follow a 
20 periodic distribution that is determined by the mass of the pair-matched 
nucleotide set, and, (v) the target nucleic acid is detected by virtue of its 
molecular weight or its sequence is determined from the mass shift of 
each extension product from its corresponding periodic reference mass. 



25 matched, the mass of each terminating base pair is unique and 

resolvable, so that the mass shifts corresponding to each terminating 
base pair are unique. The nucleotide terminators are optionally mass- 
matched or can be of distinct masses as long as distinct values of mass 
shift are obtained for each terminating base pair. 



In embodiments described above, the extending bases are pair 
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In another embodiment, the extension products are treated to 
produce blunt-ended double-stranded extension products by methods 
known to those of skill in the art, such as the use of single-strand specific 
nucleases. In an aspect of this embodiment, a plurality of target nucleic 
5 acids can be multiplexed in a single reaction by annealing each target 
nucleic acid to a primer of distinct molecular weight. The primers can be 
selected to differ in molecular weight by a value that is greater than the 
maximum mass shift, i.e., the difference in molecular weight between the 
heaviest and the lightest nucleotide terminating base pairs. Since double 

10 stranded nucleic acid can be analyzed, the effective sequence read is 
halved relative to the embodiment employing mass-matched nucleotides, 
but the number of molecules that can be multiplexed is doubled, due to 
the increase in period (the value of the mass of a base pair, rather than a 
single mass-matched nucleotide). In exemplary embodiments, about 14 

15 to about 50 sequences are multiplexed. In detection embodiments, about 
50 or more molecules can be simultaneously multiplexed since only a 
single terminating base pair is added in the extension reaction. 

In another embodiment, the chain termination reactions can each 
be carried out separately using a standard nucleotide terminator, pair- 

20 matched nucleotides, and mass-labeled primers, if modified nucleotide 
terminators which are either mass-matched or provide distinct values of 
mass shift for each terminating base pair are not available. The reactions 
can be pooled prior to detection or sequence analysis. In one 
embodiment, the mass-labeled primers can have distinct values of 

25 molecular weight that give rise to unique values of mass shift or 
positional mass difference for each terminating base. 
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Optimizing the mass spectrometric analysis of oligonucleotide mixtures 

In another method provided herein, nucleotide analogs are used to 
restrict the possible values of molecular weights that an oligonucleotide 
can possess relative to other oligonucleotides of the same length. The 
5 nucleotide analogs can be incorporated into the oligonucleotides using 
any suitable method, such as automated DNA synthesis {Oligonucleotides 
and Analogues: A Practical Approach, F. Eckstein, ed., IRL Press Oxford, 
1991) or by enzymatic replication using a polymerase and the requisite 
nucleotides and nucleotide analogs. 

10 For example, any two oligonucleotides with the same length n' 

with different base compositions can either 1) have exactly the same 
average molecular weight, or 2) have molecular weights no closer than a 
minimum value called the peak separation. In most cases, the peak 
separation will be a positive integer greater than one, but fractional values 

15 are theoretically possible. 

To illustrate an exemplary implementation of this method, the 
average molecular weight of the nucleotide analog 7-deaza-dG (328.2 
daltons) can be substituted for Mg, into equation (ii) above, which defines 
Mavg as a function of the length "n" of the oligonucleotides in bases, as 

20 follows: 

(ii) Mavg = k + n'Mc + KMt - Mc) + a(MA - Mc) + g(MG - Mc), where 
Mc, Mr, Ma, Mg are the average molecular weights of each of the four 
nucleotide bases (cytosine, thymine, adenine, guanine); c, t, a, g 
represent the number of each base present in the oligonucleotide, the sum 

25 thereof, i.e., c + t + a + g = n', the total oligonucleotide length in 
bases; and the term k is a constant representing the mass of any other 
chemical groups on the molecule, such as terminal phosphates. 
Substituting the masses of the naturally occurring bases dC, dT and dA in 
DNA (to one-tenth dalton), and of 7-deaza-dG, 

30 Mc = 289.2 
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Mt = 304.2 
Ma = 313.2 
Mg = 328.2 

and following simplification, the equation reduces to: 

Mavg = k + 289. 2n' + 15t + 24a + 39g 
Extracting the common factor from the last three terms yields 
(vi) Mavg = k + 289. 2n' + (5t + 8a + 13g) x 3 

In this example, the minimum peak separation is three daltons. It is 
not possible to identify or detect two oligonucleotides of the same length 
with different molecular weights that are closer than three daltons. 
Oligonucleotides with average masses closer than three daltons the 
oligonucleotides are detected if they are of different lengths. 

As a second example, Mt can be substituted with the molecular 
weight of a hypothetical nucleotide analog whose mass is 305.2 into 
equation (ii), yielding 

Mavg = k + 289. 2n' + 16t + 24a + 40g 
Extracting the common integer factor from the last three terms yields 

(vii) Mavg = k + 289. 2n' + (2t + 3a + 5g) x 8 
for a minimum peak separation of eight daltons. Thus, appropriate 
selection of nucleotide analogs permits construction of nucleotide sets 
that provides sufficient peak separation for adequate resolution by mass, 
such as in a time-of-flight mass spectrometer. The trade-off for a greater 
peak separation is a greater number of base compositions that have 
exactly the same mass for a given oligonucleotide length. The maximum 
number of allowed oligonucleotide masses, L, for a given nucleotide set, 
is given by 

(viii) L = n'(Mheaw - Mnght ) + 1, 
S 

where n' is the oligonucleotide length in bases, S is the peak separation, 
Mnght the mass of the lightest nucleotide in the set, Mheaw is the mass of 
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the heaviest nucleotide in the set. The number of allowed oligonucleotide 

masses scales in direct proportion to the base length and inversely with 

the peak separation, but not all possible mass values will be represented 

for a given oligonucleotide length, particularly for small n. The average 

5 density of different base compositions per allowed mass value, D, can be 

obtained by dividing equation (v) by (viii) 

D = Ntotal, 
L 

which expands into 

10 (ix) D = S(n' + l)(n f + 2(n' + 3) 

6(n'(M G - Mc) + S) 

using a typical nucleotide set with G as the heaviest base and C as the 

lightest. The density function scales in direct proportion to the peak 

separation and as a quadratic function of the oligonucleotide length in 

15 bases. In practice, the average density of base compositions per allowed 
mass value predicated by equation (ix) will be somewhat lower than the 
actual density of base compositions per observed mass value, because 
not all allowed masses will always be represented. The Examples 
describe implementation of the methods for sequencing. 

20 System and Software method for Force Mass Modulation 

Also provided are systems that automate the methods for 
determining a nucleotide sequence of a target nucleic acid or the 
detection methods provided herein using a computer programmed for 
identifying the sequence or target nucleic acid identity based upon the 

25 methods provided herein. The methods herein can be implemented, for 
example, by use of the following computer systems and using the 
following calculations, systems and methods. 

An exemplary automated testing system contains a nucleic acid 
workstation that includes an analytical instrument, such as a gel 

30 electrophoresis apparatus or a mass spectrometer or other instrument for 
determining the mass of a nucleic acid molecule in a sample, and a 
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computer capable of communicating with the analytical instrument (see, 
e.g., copending U.S. application Serial Nos. 09/285,481, 09/663,968 and 
09/836,629; see, also International PCT application No. WO 00/60361 
for exemplary automated systems). In an exemplary embodiment the 
5 computer is an IBM compatible computer system that communicates with 
the instrument using a known communication standard such as a parallel 
or serial interface. 
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For example, systems for analysis of nucleic acid samples are 
provided. The systems include a processing stations that performs a 
forced mass modulation chain extension reaction; a robotic system that 
transports the resulting products from the processing station to a mass 
5 measuring station, where the masses of the products of the reaction are 
determined; and a data analysis system, such as a computer programmed 
to identify nucleotides using forced mass modulation data, that processes 
the data from the mass measuring station to identify a nucleotide or 
plurality thereof in a sample or plurality thereof. The system can also 

10 include a control system that determines when processing at each station 
is complete and, in response, moves the sample to the next test station, 
and continuously processes samples one after another until the control 
system receives a stop instruction. 

The computer can be part of the instrument or another system 

15 component or it can beat a remote location. A computer system located 
at a site distant from the instrument can communicate wit the instrument, 
for example, through a wide area network or local area communication 
network or other suitable communication network. The system with the 
computer is programmed to automatically carry out steps of the methods 

20 herein and the requisite calculations. For embodiments that use 

mass-matched deoxyriboucleotides, a user enters the primer sequence or 
primer mass, the periodic reference mass and mass of an individual 
mass-matched deoxyonucleotide. These data can be directly entered by 
the user from a keyboard or from other computers or computer systems 

25 linked by network connection, or on removable storage medium such as a 
CD-ROM, minidisk (MD), DVD, floppy disk or other suitable storage 
medium. Next the user causes execution software that operates the 
system in which the mass spectrum of the extension products is 
generated. The Forced Mass Modulation software performs the steps of 

30 obtaining the masses of the fragments generated by the sequencing 
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reaction and measured by the analytical instrument, and determining the 
identity of a nucleotide at any base position or the positional mass 
difference. The identity of the nucleotide at each base position is 
determined by comparing the calculated Mditf[n] values to a database of 
5 previously calculated values of Mdiff for each of the chain terminating 
nucleotides. 

Mdiff [n] = Mobs[n] - MpFt[n], 

where: 

(i) MPFt[n] = (Mprimer + Mlight) + (n - 1 ) Pbase, 

10 in which n is the base position, MpR[n] is the n th periodic reference mass, 
Mpnmer is the mass of the primer, Mn g ht is the mass of the lightest 
nucleotide terminator and Pbase is the base periodicity in daltons. The 
observed masses of the sequencing reaction products are given by the 
following equation: 

15 (ii) Mobs[n] = Mpnmer + (n - 1 ) Pbase + Mterm[n], 

where n is the base position, Mobs[n] is the n th observed mass, Pbase is the 
base periodicity, and Mterm[n] is the mass of the n th terminating nucleotide 
in daltons. The positional mass differences for the sequence can be 
obtained by subtracting equation (i) from equation (ii) and evaluating at 
20 every base position n: 

where Mdiff[n] is the n th positional mass difference. This relation simplifies 

to: 

(Hi) Mdiff[n] = Mterm[n] - Mlight. 

Hence, the periodicity is determined by the mass of the mass- 
25 matched nucleotide and the shift is the difference in location of a peak 

resulting from the chain terminator. For example, in Figure 2, the lightest 
terminator is ddC, and the differential is O for C, 40 for G, 34 for A, 15 
for T. The selected mass matched nucleotide has a mass of 310 Da. The 
primer in Figure 2a has a mass of 3327 Da and the first peak would be at 
30 3600 if the first nucleotide in the extension product were C (0 shift). 
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Since the first peak is at 3640, the shift is 40 Da. Therefore the first 
nucleotide is G, corresponding to a shift from the periodicity of 310 Da 
generated by the mass-matched nucleotides. 
Detection methods 

The methods herein may be used with any method for detection of 
nucleic acids based on molecular mass known to those of skill in the art, 
particularly methods in which a primer is extended. Such methods are 
modified by extending using mass matched nucleotides and/or chain 
terminators in extension reactions. Alternatively, or additionally, 
amplification reactions may be performed using mass-matched nucleotides 
or pair-matched sets of nucleotides. These methods can be readily 
multiplexed using the methods and nucleic acid molecules provided 
herein. 

Detection methods and protocols, including those that rely on mass 
spectrometry (see, e.g., U.S. Patent No. 6,194,144; 6,225,450; 
5,691,141; 5,547,835; 6,238,871; 5,605,798; 6,043,031; 6,197,498; 
6,235,478; 6,221,601; 6,221,605; International PCT application No. WO 
99/31273, International PCT application No. WO 98/20019), can be 
modified for use with the methods herein by using mass-matched 
nucleotides for extension or pair matched duplexes for hybridization 
reactions. 

Among the methods of analysis herein are those involving the 
primer oligo base extension (PROBE) reaction with mass spectrometry for 
detection. In such reactions, the primer will be extended by mass- 
matched nucleotides. The methods herein are designed for multiplexing 
so that a plurality of different primers can be extended at different loci in 
the same reaction. The PROBE method uses a single detection primer 
followed by an oligonucleotide extension step to give products, which can 
be readily resolved by mass spectrometry, and, in particular, MALDI-TOF 
mass spectrometry. The products differ in length depending on the 
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presence or absence of a polymorphism. In this method, a detection 
primer anneals adjacent to the site of a variable nucleotide or sequence of 
nucleotides and the primer is extended using a DNA polymerase in the 
presence of one or more dideoxy NTPs and, optionally, one or more 
5 deoxy NTPs. The resulting products are resolved by MALDI-TOF mass 
spectrometry. The mass of the products as measured by MALDI-TOF 
mass spectrometry makes possible the determination of the nucleotide(s) 
present at the variable site. Use of primers containing mass-matched 
bases increases the resolving power of the reaction and permit 

10 simultaneous detection of a plurality of mutations (polymorphisms). 

These methods can be automated (see, e.g., copending U.S. 
application Serial No. 09/285,481 and published International PCT 
application No. PCT/US00/081 1 1 , which describes an automated process 
line) and performed in a system that includes a computer programmed for 

15 analysis of the mass data as described above. 

The analyses can be performed on chip based formats in which the 
target nucleic acids or primers are linked to a solid support, such as a 
silicon or silicon-coated substrate, preferably in the form of an array. 
Generally, when analyses are performed using mass spectrometry, 

20 particularly MALDI, small nanoliter volumes of sample are loaded on, such 
that the resulting spot is about, or smaller than, the size of the laser spot. 
It has been found that when this is achieved, the results from the mass 
spectrometric analysis are quantitative. The area under the signals in the 
resulting mass spectra are proportional to concentration (when normalized 

25 and corrected for background). Methods for preparing and using such 
chips are described in U.S. Patent No. 6,024,925, co-pending U.S. 
application Serial Nos. 08/786,988, 09/364,774, 09/371,150 and 
09/297,575; see, also U.S. application Serial No. PCT/US97/201 95, 
which published as WO 98/20020. Chips and kits for performing these 

30 analyses are commercially available from SEQUENOM under the 
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trademark MassARRAY. MassArray relies on the fidelity of the enzymatic 
primer extension reactions combined with the miniaturized array and 
MALDI-TOF (Matrix-Assisted Laser Desorption lonization-Time of Flight) 
mass spectrometry to deliver results rapidly. It accurately distinguishes 
5 single base changes in the size of DNA fragments associated with genetic 
variants without tags. 

The following Examples are included for illustrative purposes only 
and are not intended to limit the scope of the invention. 
TO EXAMPLE 1 

Forced Mass Modulation using Mass-Matched Deoxynucleotides 

For this implementation, a set of nucleotide analogs for the four 
bases in DNA are selected (Amersham Pharmacia Biotech) such that each 
base has exactly the same molecular weight, termed a mass-matched 

15 deoxynucleotide set. This is achieved by judiciously choosing chemical 
modifiers of the existing bases or by the using a universal base analog 
such as deoxyinosine, which can form base pairs with more than one of 
the natural bases. For this example, the mass of each deoxynucleotide 
("dN") in the mass-matched set has the arbitrarily selected value of 310 

20 daltons, but any other value suffices. The sequencing reaction is 
performed as follows: 1) a primer is annealed to the target to be 
sequenced; 2) the resulting structure is subjected to a extension reaction 
using a suitable polymerase in the presence of the mass-matched 
nucleotide set and the four standard dideoxynucleotide terminators. The 

25 products and molecular masses of such a reaction are shown with a 

simulated mass spectrum in Figure 2a. The base periodicity is the mass 
of dN, or 310 daltons. The identity of a nucleotide at any base position is 
given by the positional mass difference, defined as the distance in daltons 
between the observed peak and the nearest periodic reference mass, 

30 which occurs every 310 daltons. In this example, the first periodic 
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reference mass is defined as the (primer mass + ddC), or (3327 + 273) 
= 3600 daltons. The second periodic reference mass would be 3600 
plus the base periodicity or (3600 + 310) = 3910, and so on. 
Expressed in terms of the base position n: 

(j) MpR[n] = (Mpr.mer + Mlight) + (n - 1 ) Pbase, 

where n is the base position, MpR[n] is the n th periodic reference mass, 
Mpnmer is the mass of the primer, Mnght is the mass of the lightest 
nucleotide terminator and Pbase is the base periodicity in daltons. The 
observed masses of the sequencing reaction products are given by the 
following equation: 

(ii) Mobs[n] = Mpnmer + (n - 1) rbase + Mterm[n], 

where n is the base position, IVUs[n] is the n th observed mass, Pbase is the 
base periodicity, and IvWn] is the mass of the n th terminating nucleotide 
in daltons. The positional mass differences for the sequence can be 
obtained by subtracting equation (i) from equation (ii) and evaluating at 
every base position n: 

Mdiff[n] = Mobs[n] - MpR[n], 
where Mdi«[n] is the n th positional mass difference. This relation simplifies 
to: 

(ii'l) Mdiff[n] = Mterm[n] - Mlight. 

Inspection of equation (iii) reveals that Mdiff can only take on four distinct 
values, each corresponding to a different nucleotide terminator: 
Md,ff["ddC"] = (273.2 - 273.2) = 0 
Md,ff["ddT"] = (288.2 - 273.2) = 15 
Mdif f ["ddA"] = (297.2 - 273.2) = 24 
Md,ff["ddG"] = (313.2 - 273.2) = 40. 
Hence, the identity of the nucleotide at every base position in the target 
sequence can be determined by comparing each calculated positional 
mass difference with the values in the table above. Since the values that 
Mdiff can assume depend only on the choice of nucleotide terminators 
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used in the sequencing reaction, it is possible to tailor the positional mass 
differences so that they are resolvable for any particular mass 
spectrometer. For example, replacing the terminator ddT with its analog 
5-bromo-dideoxyuridine (353.1 daltons) yields a positional mass 
5 difference of (353.1 - 273.2) = 79.9 Da for termination at T positions in 
the target sequence. This type of nucleotide substitution can be 
particularly valuable for lower-resolution mass spectrometers, as it 
possible to maintain the sequence read accuracy without requiring any 
additional mass spectra. 

10 Further inspection of equation (iii) reveals that each observed mass 

value can be at most 40 daltons heavier than the nearest periodic 
reference mass. This limit is termed the maximum mass shift and is 
defined as the mass difference between the heaviest nucleotide 
terminator and the lightest. Resolving a second target sequence by 

15 Forced Mass Modulation with the standard dideoxy terminators is possible 
in a single spectrum so long as the primer for the second sequence is at 
least 40 daltons heavier (the maximum mass shift) than the primer for the 
first sequence, thus insuring that the peaks for each sequence never 
overlap in mass. 

20 In practice, it is recommended most mass spectrometric formats 

that the second primer is at least about 60 daltons heavier than the first 
primer, as each observed peak will have a finite width. Figure 2b shows 
a target second sequence resolved on the same mass spectrum shown in 
Figure 2a, using a primer heavier by 77 daltons. The peaks corresponding 

25 to the reaction products from the first target sequence can fall within the 
shaded regions of the spectrum in Figure 2b, which can never intersect 
peaks from the second target sequence. Unambiguous resolution of both 
sequences is possible in this arrangement because each peak can be 
uniquely assigned to a nucleotide, a base position, and a target sequence. 

30 This method is designated Mass Spectrum Division Multiplexing herein, 
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and it is implemented using mass-staggered primers. Figure 3 shows four 
different sequences resolved in a single spectrum using a set of 
mass-staggered primers that are separated in mass by integer multiples of 
77 daltons (77, 154, and 231 daltons). 
5 The theoretical upper limit on the number of sequences that can be 

multiplexed in a single mass spectrum is given by the following equation: 

(iv) L' = Pbase 
Smax, 

where L' is the upper limit, Pbase is the base periodicity, and Smax is the 

10 maximum mass shift in daltons. For the nucleotide set and terminators 
used in this example, L = {310 / 40) = 7.75, or approximately seven. 
Increasing the number of sequences that can be multiplexed in a single 
spectrum, can be achieved by implementing one or both of an increase in 
the base periodicity, and a reduction of the maximum mass shift. The 

15 base periodicity can be increased by choosing a mass-matched nucleotide 
set that has a higher molecular weight for dN. It is simpler to lower the 
maximum mass shift by careful use of the nucleotide terminators and their 
analogs. For example, if the sequencing reactions were performed using 
only the terminators ddC, ddT, and ddA, then the maximum mass shift 

20 becomes (mass of ddA - mass of ddC) = (297 - 273) = 24 Da. In this 
case the upper limit on the number of sequences that can be multiplexed 
is L = (310 / 24) = 12.92, or approximately twelve. In situations where 
complete sequence information is not required, such as diagnostic 
sequencing, a great reduction in the number of required spectra can be 

25 realized by using fewer than four nucleotide terminators. If the 
sequencing reaction is performed using only a single nucleotide 
terminator, the maximum mass shift becomes identically zero, and the 
number of sequences that can be multiplexed in a single spectrum is 
limited only by the absolute resolution of the mass spectrometer in 

30 question. If a given mass spectrometer has an absolute resolution of 12 
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Da in the mass range of the sequencing reaction products, then the 
maximum number of sequences that can be multiplexed is given by L = 
(310 / 12) = 25.83, or approximately twenty-five. 

EXAMPLE 2 

5 Forced Mass Modulation using Pair-Matched Deoxynucleotides 

Implementation of Forced Mass Modulation using pair-matched 
nucleotides is shown in Figure 4. The basic requirement for this method 
is that the sequencing reaction products can be analyzed as 
double-stranded structures. Briefly, the steps in the reaction are as 

10 follows: 1) A partially duplex hairpin primer with a 3' overhang and a 5' 
phosphate group is annealed and ligated to the single stranded target 
sequence. 2) The resulting partially duplex structure is subjected to a 
sequencing reaction using the pair-matched nucleotide set described 
above along with the set of mass-matched terminators {ddM). 3) The 

15 products from the sequencing reaction are exposed to a strict single 
strand-specific nuclease that results in the production of blunt-ended 
hairpin structures ready for analysis by mass spectrometry. Figure 5 
shows the products and molecular masses of the nuclease digestion along 
with a simulated mass spectrum. 

20 Because the reaction products are double-stranded, they are forced 

to assume a quasi-periodic distribution with a base periodicity of 617.4 
daltons. The shaded regions on the spectrum shown in Figure 5 indicate 
the allowed mass ranges that can be occupied by the reaction products. 
The first periodic reference mass is at 10360 Da, which is the mass of 

25 the fully duplex hairpin primer plus a ddM:dC base pair. Expressing the 
periodic reference masses in terms of the base position n yields: 

(X) MpR[n] = (Mduplex + Mlight + MddlVl) + (n - 1 ) X Pbase 

Where Mppt[n] is the n th periodic reference mass, Mduplex is the mass of the 
fully duplex primer, Mhght is the mass of the lightest deoxynucleotide in the 
30 target, Pbase is the base periodicity, and MddM is the mass of ddM in 

-49- 



25491-2408 



daltons. The observed masses of the sequencing reaction products are 
given by the following equation: 

(Xl") Mobs[n] = Mduplex 4- MddM + (n - 1) X Pbase + Mtarg[n], 

where n is the base position, Mobs[n] is the n th observed mass, and Mtar g [n] 
5 is the mass of the n th nucleotide in the target sequence past the priming 
site in the 3' -> 5' direction. 

In contrast to the mass-matched nucleotide set implementation that 
provides the sequence complementary to the template strand read in the 
5' -> 3' direction, the pair-matched nucleotide set implementation 

10 described herein directly reads the template strand in the 3' -> 5' 

direction. The positional mass differences for this implementation are the 
same as those in Example 1 , except that the mass difference 
corresponding to a termination on dG is 39 as opposed to 40 daltons, 
because 7-deaza-dG is exactly one dalton lighter than dG. Since double 

15 stranded DNA can be analyzed for this method to work, the effective 
sequence read length is halved, although the number of sequences that 
can be multiplexed is doubled, due to the increase in the base periodicity. 

As a demonstration of Forced Mass Modulation implemented 
without using mass-matched terminators, the positional mass differences 

20 for the above example using the following set of nucleotide terminators is 
calculated as follows: 

Terminator Nucleotide Ana lop Mass Base Pairing Mass of Base 

Pair 

T 5-Bromo~dideoxyuridine 353. 1 5-Br-ddUidA 
25 666.3 

C 5-Methyl-dideoxycytidine 287.2 5-Me-ddC: 7-deaza-dG 

615.4 

A Dideoxyadenosine 297.2 ddA: dT 

601.4 

30 G Dideoxyinosine 298.2 ddl: dC 587.4 
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The positional mass difference at every base position is given by: 

(Xii) Mdiff[n] = Mpairfn] - Mlightest, 

where Mdiffln] is the n th positional mass difference, M P air[n] is the mass of 
n th terminating base pair, and Mugmest is the mass of the lightest 
5 terminating base pair in daltons. Substituting in the values from the table 
above yields: 

Mdiff["G M ] = (587.4- 587-4) = 0 
Mdiff["A"] = (601.4- 587-4) = 14 
Mdiff[ ,r C n ] = (615.4- 587.4) - 28 
10 Mdrff["T"] = (666.3 - 587-4) = 78.9 

Since each terminating base pair has a unique positional mass difference, 
the base sequence can be determined unambiguously. The maximum 
mass shift in this case is 78.9 daltons. When choosing a set of 
terminating nucleotides it is important to select the set such that the 
15 positional mass difference for each base termination is distinct and 
resolvable by mass. 

If modified nucleotide terminators are not used, it is still possible 
to implement Forced Mass Modulation by carrying out each of the four 
termination reactions separately using mass-labeled primers rather than 
20 modified terminators, combining all reaction products, and then obtaining 
a mass spectrum. In order to produce the same positional mass 
differences as shown in Example 1 , using a set of pair-matched 
nucleotides and the standard dideoxy terminators, the following primer 
mass shifts are required: 
25 Termination Reaction Primer Mass 



C "reference " primer 

T reference primer + 15 Da 

A reference primer + 24 Da 

G reference primer + 39 Da 
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This method is essentially equivalent to multiplexing four single-nucleotide 
sequencing reactions in the same spectrum, except that all the 
sequencing products originate from the same priming site but terminate 
on different nucleotides. 



Forced Mass Modulation in the Detection and Scoring of Single Nucleotide 
Polymorphisms 

Forced Mass Modulation can be used to simplify the analysis of 
closely related sequence variants, as is required in the detection and 

10 scoring of single nucleotide polymorphisms. Figure 6 shows three 
sequence variants that differ from each other only at a single base 
position sequenced by a conventional Sanger reaction. The mass 
distribution of the reaction products is so complex that it can be 
uninterpretable, even if the base sequences of the variants are known a 

1 5 priori. 

Figure 7 shows the same three variants sequenced by Forced Mass 
Modulation using mass-matched deoxynucleotides (dN = 310 Da) and the 
standard dideoxy terminators. The positions and identities of the 
single-nucleotide changes are immediately apparent from the mass 
20 spectrum. Since the masses of the sequencing reaction products are 

constrained to fall within the shaded regions of the spectrum in Figure 7b, 
it is possible to multiplex other sequences on the same spectrum. 



Base Composition Density Distributions for the Total Set of possible 7- 
25 base Oligonucleotides 

For this implementation, three sets of 7-base oligonucleotides 

comprising all possible base compositions for a 7-base oligonucleotide can 

be obtained; the first set comprising the four natural bases (dA, dG, dC 

and dT), the second set comprising three of the natural bases (dA, dC 

30 and dT) and the nucleotide analog 7-deaza-deoxyguanosine (7-deaza-dG) 



5 



EXAMPLE 3 



EXAMPLE 4 
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substituted for dG, and the third set comprising three of the natural bases 
(dA, dG and dC) and the nucleotide analog deutero-deoxythymine 
(deutero-dT) substituted for dT. Figure 8 shows the actual base 
composition density distributions for the total set of possible 7-base 
5 oligonucleotides using the three different nucleotide sets. Note that for 
the set of naturally occurring bases (Figure 8a), nearly every base 
composition has its own distinct mass value, but most of these mass 
values are spaced only one dalton from each other. Increasing the peak 
separation to three daltons by substitution of dG with 7-deaza-dG (Figure 

10 8b) markedly increases the average number of base compositions per 
observed mass, particularly for those masses in the center of the range, 
but any two oligonucleotides of the same length with different molecular 
weights will have to be separated by at least three daltons. Similarly, 
substitution of dT with deutero-dT (Figure 8c) gives a minimum peak 

15 separation between oligonucleotides having the same length but different 
molecular weights of eight daltons. The trade-off for a greater peak 
separation is a greater number of oligonucleotides that have exactly the 
same mass for a given oligonucleotide length. 

20 Since modifications will be apparent to those of skill in this art, it is 

intended that this invention be limited only by the scope of the appended 
claims. 
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